Negative differential conductance and magnetoresistance oscillations due to spin 
accumulation in ferromagnetic double-island devices 
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Spin-dependent electronic transport in magnetic double-island devices is considered theoretically 
in the sequential tunneling regime. Electric current and tunnel magnetoresistance are analyzed as a 
function of the bias voltage and spin relaxation time in the islands. It is shown that the interplay of 
spin accumulation on the islands and charging effects leads to periodic modification of the differential 
conductance and tunnel magnetoresistance. For a sufficiently long spin relaxation time, the mod- 
ulations are associated with periodic oscillations of the sign of both the tunnel magnetoresistance 
and differential conductance. 
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Introduction - Spin-dependent transport in single- 
electron devices is attracting currently much attention 
from both fundamental and application points of view 
QHHHIEHIzlH It has been shown that the interplay 
of single-electron charging effects and spin dependence 
of tunneling processes (caused by ferromagnetism of ex- 
ternal and/or central electrodes) in single-electron tran- 
sistors gives rise to novel interesting phenomena, such 
as periodic modulation of the tunnel magnetoresistance 
(TMR) with increasing either bias or gate voltages, quasi- 
oscillatory behavior of the spin accumulation on the cen- 
tral electrode (referred to as an island in the following), 
enhancement of TMR in the Coulomb blockade regime, 
and others. Some of the theoretically predicted phenom- 
ena have already been confirmed experimentally. 

Recently, transport characteristics of granular mag- 
netic nanobridges connecting two external magnetic 
and/or nonmagnetic electrodes were investigated experi- 
mentally H El Ell. 

The nanobridges consisted of mag- 
netic grains distributed in a nonconducting matrix. The 
corresponding current- voltage (TV) curves revealed char- 
acteristic steps due to discrete charging of the grains with 
single electrons (Coulomb steps). Apart from this, oscil- 
lations of the TMR effect (associated with the transition 
from one magnetic configuration to the other) with the 
bias voltage have also been observed [Tj. Moreover, two 
additional features of the transport characteristics were 
found, whose physical mechanism and origin needs fur- 
ther considerations. First, the differential conductance 
was found to change sign periodically with increasing bias 
voltage 0, [ljj • Second, the TMR was shown to oscillate 
between positive and negative values with increasing bias 
voltage. Physical origin of these oscillations, however, re- 
mained unclear, although the role of spin accumulation 
was invoked following the results obtained for a single- 
island system [2|. 

For small lateral and vertical dimensions of the 
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nanobridges, electronic transport between the external 
leads occurs as a result of consecutive electron jumps via 
metallic grains, and the corresponding electronic paths 
can include either one, two, or more grains. Transport 
through a single grain and the associated TMR have al- 
ready been analyzed 2], but the results cannot describe 
the above mentioned data obtained on magnetic multi- 
grain nanobridges. However, the experimental setup de- 
scribed above can be modelled by assuming only two 
grains involved in the transport processes |9(. Although 
the Coulomb steps have been accounted for by the ex- 
isting theoretical models based on discrete charging with 
single electrons, the origin of negative differential con- 
ductance (NDC) as well as the oscillations in TMR with 
alternating sign remains unclear. Therefore, this problem 
is considered in this paper in more details. 

Following the above discussion, we consider transport 
through a double-island device consisting of two metal- 
lic grains attached to two external leads. The two is- 
lands are separated from each other and also from the 
external leads by tunnel barriers. In a general case all 
the four metallic components of the device (two external 
electrodes and two islands) can be ferromagnetic Q, || , 
but in the following we will consider the situation with 
one external lead being nonmagnetic and the other three 
components being ferromagnetic. We also restrict our 
considerations to collinear [parallel (P) and antiparallcl 
(AP)] alignments of the magnetic moments, as indicated 
in the inset of Fig.^d). 

The key role in our analysis is played by the nonequi- 
librium spin accumulation, which occurs for a sufficiently 
long spin relaxation time on the islands. Strictly speak- 
ing, the accumulation may take place when the spin re- 
laxation time is significantly longer than the time be- 
tween successive tunneling processes. Assuming this is 
the case, we show that the interplay of charging effects 
and spin accumulation gives rise to both effects described 
above, i.e., to the NDC and periodic change of the sign 
of TMR with increasing bias voltage. This behavior of 
transport characteristics accounts for the experimental 
observations. We also show that the presence of NDC 
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depends on the magnetic configuration of the device and 
by switching from one configuration to the other one 
may change sign of the differential conductance (eg., by 
an external magnetic field). This behavior may be of 
some interest for possible future applications in spintron- 
ics / magnetoelectronics devices. 

Model - The system adapted to model the phenomena 
discussed above consists of the left and right leads and 
two central islands, separated from each other and from 
the leads by tunnel barriers, as shown schematically in 
the inset of Fig. ^d). The bias voltage is applied to 
the external leads - Vl to the left one and Vr to the 
right one. In a general case one can attach capacitively 
gate voltages to both islands. However, to model the 
experimental results described above, we neglect the gate 
voltages in the following. We assume that the size of each 
island is sufficiently small to have the charging energies 
of the islands significantly larger than the thermal energy 
fcfjT, but still large enough to neglect size quantization. 
Apart from this, we assume that the barrier resistances 
are much larger than the quantum resistance, R r ;§> h/e 2 
(r = L, M, R). The system is then in a well defined charge 
state described by ri\ and n-i excess electrons on the first 
(left) and second (right) islands, respectively, and the 
electrostatic energy of the system is given by the formula 

E3,m 



E(ni,n 2 ) = E Cl ni 



Ec 2 n 2 



CrVr 
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+2E Cm I ni — J I 112 — J , (1) 

where Cl(r) is the capacitance of the left (right) junc- 
tion, Ec 1 and Eq 2 denote the charging energies of the 
two islands, £ Cl( c 2 ) = e 2 /(2C 1(2) )[l - C^/ (C x C 2 )]-\ 
and Eq m is the energy of electrostatic coupling between 
the islands, E Cm = e 2 / (2C M )(C 1 C 2 /C 2 l - l)" 1 , with 
Ci(2) being the total capacitance of the first (second) 



island, C*i(2) = C] 



^l(r) + Cm i and Cm denoting the ca- 
pacitance of the middle junction (the one between the 
islands). This allows us to employ the quasi-classical the- 
ory based on the master e quat ion and the Fermi golden 
rule for the tunneling rates [lj] . Such approach describes 
well transport in the sequential tunneling regime and cor- 
responds to taking into account only the first-order tun- 
neling processes, which are exponentially suppressed in 
the Coulomb blockade regime but give the dominant con- 
tribution to electric current when the applied voltage ex- 
ceeds a certain threshold. Moreover, we take into account 
only non-spin-flip (spin-conserving) tunneling processes 
through the barriers. 

When a bias voltage V = Vl — Vr is applied to the sys- 
tem and spin relaxation time on the islands is sufficiently 
long, a nonequilibrium magnetic moment can accumulate 
on each island. Let us denote the corresponding shift of 
the Fermi level for electrons with spin-cx on the j-th island 
by A-Epj . When in the initial state there were n\ and n 2 
excess electrons on the islands, the spin-dependent tun- 



neling rate from the left electrode to the first island is 
then given by the formula 



1^(711,712) 
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AE^( ni ,n 2 ) 



e 2 Rl cxp [AEl 1 (n 1 ,n 2 )/k B T] - 1 
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where i?£ is the spin-dependent resistance of the left bar- 
rier and Ai?Li {n\ , 772) describes the change in the electro- 
static energy of the system caused by the respective tun- 
neling event, AS£ 1 (ni, n 2 ) = E(ni + 1, 7i 2 ) - E(ni, n 2 ) + 
eVL + AE^ 1 - (According to our notation the electron 
charge is — e with e > 0.) 

For the following discussion it is convenient to distin- 
guish between spin projection on the global quantiza- 
tion axis (cr =], I) an( i on t ne local quantization axis 
(a = +,— , with a — + and a = — corresponding to 
spin-majority and spin-minority electrons, respectively). 
The spin asymmetry in tunneling processes follows from 
spin-dependent density of states and spin-dependent tun- 
neling matrix elements. Let us define the parameters 
P r = D+/D- for the leads (r = L, R) and fa = D+fD^ 

for the islands (j = 1,2), where Dr an d D^' are the 
appropriate densities of states for spin-majority (spin- 
minority) electrons in the leads (r = L, R) and islands 
(j = 1,2), respectively. Thus, the Fermi level shifts due 
to spin accumulation on the islands obey the condition 
AEp-/AEp- — ~f3j for the first (j = 1) and second 
(j — 2) islands. 

The spin asymmetry of the barrier resistances can be 
described by the parameters a r = R}. / Rj. for r = L, M, R. 
Assuming a constant (independent of energy) density of 
states and constant (independent of energy and spin) ma- 
trix elements, one can write a r = D^D^-/ D],d\- for the 
left (r = L, j = 1) and right (r = R, j = 2) barri- 



ers, and a r = D^D^/ D^D^ for the central (r = M) 
barrier. When the local and global spin quantization 
axes coincide (parallel configuration), then one can write 
ay = l/(p r f3j) for the left (r = L, j = 1) and right 
(r = R, j = 2) barriers, and a r = l/{(3\(3 2 ) for the cen- 
tral (r = M) barrier. The above formulas are also appli- 
cable to the situation with magnetic moment of a lead or 
an island reversed (antiparallel configuration), but with 
the corresponding (3 replaced by 1//3. 

The probability P(ni,7i2) that the system is in a 
charge state (771,772) can be determined in a recursive 
way from the appropriate steady-state master equation 
|13|. The electric current flowing through the left junc- 
tion is then given by 

^L = -eJ^ K 1 (m,n 2 ) -r^ L (ni,n 2 )]P(ni,n 2 ) .(3) 

The associated shifts of the Fermi level can be calculated 
in a self-consistent way from the relations 0] 



M(R) 1 L(M) 



- P n(2)°ii(2) 

T sf,l(2) 



A^ 1(2) =0, (4) 



where flij is the volume of the island j, r s f.j denotes the 
spin-flip relaxation time in the j-th island, and 1° is the 
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lei and antiparallel configurations are denned. This sys- 
tem geometry corresponds to the situation studied ex- 
perimentally in Ref. @. The transport characteristics 
have been obtained for the limit of long spin relaxation 
time. The current-induced shifts of the Fermi level (due 
to spin accumulation) for the spin-majority electrons in 
both islands are shown in Fig.^a) for the P and AP con- 
figurations and for positive bias voltage (electrons flow 
from right to left). A nonzero spin accumulation occurs 
in both magnetic configurations. In the P configuration 
the shift of the Fermi level for spin-majority electrons is 
negative for both islands, whereas in the antiparallel con- 
figuration it is positive for the first island and negative 
for the second one. Such a behavior can be accounted 
for by taking into account spin asymmetries of the tun- 
neling processes through all the three barriers, similarly 
as it was done in the case of a double-barrier system 0]. 

The Coulomb steps in the P and AP configurations 
are significantly different [see Fig. QJb)]. There are two 
reasons of this difference. First, the overall resistance 
of a given barrier depends on the relative orientation of 
the magnetic moments of adjacent ferromagnetic compo- 
nents of the device. Second, the spin accumulations in 
the P and AP configurations are also different. The latter 
fact is of particular importance for the present analysis. 

The differential conductance corresponding to the I—V 
curves from Fig.^b) is shown in Fig.^c). In both P and 
AP configurations the differential conductance changes 
sign periodically with increasing bias. However, the bias 
voltage range of NDC for the P configuration is different 
from that for the AP one. The corresponding phase dif- 
ference varies with the bias voltage. Furthermore, NDC 
is more pronounced in the AP configuration and its abso- 
lute magnitude increases with increasing bias voltage, as 
shown in Fig. ^c). This behavior of NDC is consistent 
with experimental observation reported in Ref. 0- 

The TMR effect associated with the transition from 
antiparallel to parallel configurations is visualized in 



FIG. 1: Shifts of the Fermi levels for spin-majority electrons 
(a) and currents (b) in the P and AP configurations; differ- 
ential conductance in the AP configuration (c); and TMR 
(d) as a function of the bias voltage for r B f t i = 7-^2 — > 00. 
The other parameters are: T = 140 K, Cl = 0.45 aF, 
Cm = 0.2 aF, C R = 0.35 aF, fa = = Pr = 0.2, and 
/3l = 1. The total barrier resistances in the P configuration 
are: -Rl/3500 = -Rm = Rr = 1 Mf2. In the AP configuration, 
R$ = Rr* = (Rr,iRr,i) 1/2 , for r = M,R. The parameters 
correspond to the experimental ones taken from Ref. |J 



current flowing through the barrier r (r = L, M, R) in the 
spin channel a. 

Results and discussion - In Figure ^ we show trans- 
port characteristics of a device whose two islands as well 
as the right electrode are made of the same ferromag- 
netic metal, whereas the left electrode is nonmagnetic, 
as shown in the inset of Fig.^d), where also both paral- 
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FIG. 2: (color online) Differential conductance in the antipar- 
allel configuration as a function of the bias voltage calculated 
for different values of the spin relaxation time r s f,i = r s f,2 = 
T a f and for D^fin = _D^f2i2 = 1000/eV. The other parame- 
ters are the same as in Fig. Q 



4 



0.20 




0.15 


A /"\ T si = 10_9s - 

/ \ < ■ limit of short , _-7 

f.\ 1 1 , „ u 1 =10 s 

Ml ■ . spin relaxation st 


0.10 
0.05 
> u.uu 


M ' ■■' / X -10 S 

/ 'i M — -t%io"s J 
/ '• V-yV'- \sJ's-\ C-\ -"A 

* ' -A p- / J^- jS^HJH \^4c j 


\ ; 1 i-i > r i f xJ 

->.:J H W v 


H -0.05 


-0.10 


\ "i n 
- \ i 


-0.15 


.' limit of long 
spin relaxation 


-0.20 





V(V) 

FIG. 3: (color online) The TMR effect as a function of the 
bias voltage calculated for different values of spin relaxation 
time T B f,i = r s f,2 = T s f and for Dj+fin = Djt,f2i2 = 1000/eV. 
The other parameters are the same as in Fig. 

Fig. _Hd) . The TMR is periodically modulated with in- 
creasing bias voltage. The initial phase of the oscillations 
depends on the system parameters and can change from 
negative to positive. Moreover, the modulations are as- 
sociated with a periodic change of TMR between positive 
and negative values. 

Oscillations of the sign of TMR and differential con- 
ductance result from spin accumulations in both islands 
and are absent in the limit of fast spin relaxation (no 
spin accumulation) on the islands. The results presented 
in Fig. H have been calculated for the long spin relax- 
ation limit. In such a limit some spin accumulation may 
occur even for a very small current flowing through the 
system, giving rise to NDC and oscillations of the TMR 
sign for small bias voltages. However, both effects dis- 
appear when the spin relaxation time is shorter than the 
time between successive tunneling events. Thus, for a 
finite relaxation time, one may expect no NDC and no 
TMR oscillations for small voltages and an onset of these 



effects at larger voltages. This is because at some voltage 
there is a crossover from the fast to slow spin relaxation 
limits. In fact such a behavior is consistent with experi- 
mental data presented in Ref. @- 

The disappearance of NDC with decreasing spin relax- 
ation time T a t is shown explicitly in Fig. [21 where the bias 
dependence of differential conductance is presented for 
different values of the spin relaxation time. This figure 
clearly shows that NDC disappears when spin relaxation 
time decreases, in agreement with the above discussion. 
In the limit of fast spin relaxation time, the differential 
conductance is positive, although its periodic modulation 
still remains. 

Similarly, periodic oscillations of the sign of TMR also 
disappear with decreasing spin relaxation time r s f. This 
behavior is shown in Fig. [3J where the bias dependence 
of TMR is shown for several values of r s f. First, the 
transitions to negative TMR disappear with decreasing 
r s f. The TMR becomes then positive, although some 
periodic modulations survive. Second, the phase of the 
modulations shifts by about tt when the spin relaxation 
varies from fast to slow limits. 

In conclusion, we have analyzed transport through 
double-island device and shown that the negative dif- 
ferential conductance measured experimentally is due to 
nonequilibrium spin accumulation in the islands. Fur- 
thermore, spin accumulation may also lead to oscillations 
in TMR between negative and positive values. The effect 
of NDC occurs in both configurations and transition from 
one configuration to the other may result in transition 
from positive to negative differential conductance, which 
may be of some importance for applications in spintronics 
devices. 
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